Quantifying hierarchical mixture quality in polymer composite materials: structure 

and inhomogeneity in multiple scales 
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Mixture quality plays a crucial role in the physical properties of multi-component immiscible 
polymer mixtures including nanocomposites and polymer blends. Such complex mixtures are often 
characterized by hierarchical internal structures, which have not been accounted for by conventional 
mixture quantifications. We propose a way to characterize the mixture quality of complex mixtures 
with hierarchical structures. Starting from a concentration field, which can be typically obtained 
from TEM/SEM images, the distribution of the coarse-grained concentration is analyzed to obtain 
the scale-dependent inhomogeneity of a mixture. The hierarchical nature of a mixture is character- 
ized by multiple characteristic scales of the scale-dependent inhomogeneity. We demonstrate how 
the proposed method works to characterize sizes and distributions in different dispersions. This 
method is generally applicable to various complex mixtures. 

Keywords: Mixing, Polymer composites, Filler dispersion, Image analysis, Multiple scales of segregation, 
Multi-scale measure 



I. INTRODUCTION 

Multi-component polymer materials including 
nanocomposites, fiber-reinforced plastics, and poly- 
mer blends typically have inhomogeneous internal 
structures on a mesoscopic level. The properties of 
such mixture materials are strongly dependent on the 
distribution of their internal structures, namely, the 
mixture quality Therefore, there is a common 

need to develop a quantification measure to describe the 
mixture quality of multi-component polymer materials. 
Such a technique would also be required to evaluate 
mixing processes in various mixing devices in the field 
of chemical engineering Though various methods 

for the analysis of mixture quality have been proposed, 
these approaches are sometimes system-specific and 
fail to directly assess hierarchical internal structures. 
To compare different mixture systems (different filler 
loading, different mixing conditions, etc.), a simple 
quantification method based on non-system-specific 
criteria is required. 

In general, internal structures in a complex mixture 
have their own characteristic scales. Therefore, the si- 
multaneous quantification of size and distribution is nec- 
essary. However, deviation from uniformity, which is a 
commonly used method, often misses characteristic scales 
of internal inhomogeneity because uniform distribution 
itself is intrinsically defined without any scale. From this 
observation, the identification of characteristic scales is a 
fundamental problem when discussing the mixture qual- 
ity of complex mixtures. 

Many prior works regarding the quantification of mix- 
ture quality primarily focused on the non-uniformity of 
distribution, regardless of the sizes or structures of the 
dispersed phase. Typically for a dispersion of spherical 



particles, the center points of the particles are first identi- 
fied, then the distribution of the center points is analyzed 
|15l - il7| . One direct approach to measure the distribu- 
tion inhomogeneity is to evaluate the deviation from the 
uniform distribution based on fluctuation of the point 
density [H, EE fl8rl2^ |. Another approach is to evalu- 
ate a certain cost function of the inter-point distance, 
which is the minimum when the points are uniformly 
distributed [17] . These techniques are useful for systems 
where the ideal distribution is homogeneous but not for 
systems with internal structures. The length scales of 
internal structures cannot assessed by these approaches. 

Instead of measuring deviation from uniformity, a de- 
gree of clustering was used to characterize the inhomo- 
geneity of the distribution. The spatial correlation func- 
tion of density of the dispersed phase is the simplest 
way to define correlation length, which is an average 
size of the cluster d, [23|. Alternative methods to de- 
fine size were developed based on either a concentration 
gradient [Hj], or the inter-point distance [ll|. The size 
of the matrix phase not containing the dispersed phase 
can also be used to characterize the clustering tendency, 
which is closely related to the mechanical reinforcement 
of nanocomposites [24j. Another method is based on 
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the mathematical morphology, which considers the vol- 
ume fraction of isotropic virtual dilation of the dispersed 
phase. The dependence of the volume fraction on dilation 
differs depending on the degree of clustering [25| . These 
approaches focus on one representative length scale but 
not on the multiple characteristic scales associated with 
hierarchical internal structures. 

Internal structures within a complex mixture are often 
hierarchical in nature. Thus, different scales associated 
with hierarchical structures should be involved in charac- 
terizing the mixture quality of a complex mixture. The 
only other known study that has assessed inhomogene- 
ity at different scales involves the Mix-Norm method [26| . 
However, the Mix-Norm method has focuses on the defin- 
ing non-uniformity of a system in which the ideal mixture 
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state is uniform for all scales and thus does not apply to 
the identification of hierarchical characteristic scales. 

In this paper, we propose a scale-dependent measure 
for inhomogeneity of complex mixtures that is capable of 
identifying multiple characteristic scales of inhomogene- 
ity. The degrees of inhomogeneity are defined at different 
scales through the coarse-graining of the concentration 
field of a dispersed phase. Multiple characteristic scales 
are identified with the scale dependence of inhomogene- 
ity. We demonstrate how the scale-dependent measure 
works by applying the technique to synthetic dispersions 
with different distributions. This technique can aid in 
comparing different complex mixtures and can provide 
a fundamental understanding of the mixture quality of 
systems with different internal structures. 



II. METHOD OF SCALE-DEPENDENT 
MOMENT FOR MIXTURE QUALITY 

We start with a concentration field, c(x), of a com- 
ponent in a mixture. The inhomogeneity of the mixture 
can be characterized based on the fluctuation in the con- 
centration field, which contains information on a micro- 
scopic level. The resolution of c(x) is determined by the 
resolution of the data obtained. 

Because the definition of inhomogeneity should not de- 
pend on the total content of the component, we consider 
the normalized concentration field, namely a probability 
density, as 



c(x) 



In A v<y) ' 



(1) 



where Q denotes the whole domain of the mixture. The 
induced probability density fi(x) satisfies J n dxfi(x) = 1 
irrespective of the total content of the dispersed sub- 
stance J Q dxc(x). 

To quantify the fluctuation of a specific scale, the 
coarse-grained field is defined as 
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where B r {x) denotes the domain of a linear size r 
around a location x. Coarse- graining is a familiar con- 
cept in statistical physics and has been especially use- 
ful in characterizing fractal structure in fluid turbulence 
[27l - l30l |. Coarse-graining has also been applied to char- 
acterize time-series with long-range correlation [3lj . In 
these applications, fractal structures were characterized 
by power-law scaling, and the exponents of this analy- 
sis were a primary concern. Here, we are interested in 
applying the coarse-graining concept to quantification of 
mixture quality. 

The fluctuation at a resolution scale r is characterized 
by the statistical moment, 



(3) 



{i4) = / dxi4(x), q e (-00,00), 



where the parameter q represents the magnitude of fluc- 
tuation: a large positive q indicates a large fluctuation 
and a negative q indicates a small fluctuation. 
If a mixture is homogeneous at a scale r, 



= for q^P, 



(4) 



holds. For instance, in the limit of macroscopic scale, 
r — ► L, where L denotes the system size, and the coarse- 
grained concentration becomes \xl = J Bl dxfi(x) « 
Jq dx[i(x) by definition; therefore, the normalized mo- 
ment is always close to unity at the system size, 
~ This property simply represents the 
macroscopic homogeneity in the general situation. 

Furthermore, if there is no coherent structure of scale 
r, the relationship 



( M <?) ~ Ar D \ 



(5) 



holds with a coefficient A, where D is the spatial dimen- 
sion. Based on the properties (j4j and ((5|) , to characterize 
the inhomogeneity of a mixture, we define the inhomo- 
geneity function at scale r by a normalized moment, 



F q {r) = 



ML 



1/(9-1) 



(6) 



To understand the physical implication of F q (r), we 
consider a mixture with a correlation length I. In this 
case, the inhomogeneity function behaves as 



cst.(> 1) r < I, 



F q (r) { > I 



1 



r > I, 



(7) 



The fluctuation in a smaller scale than r is accounted 
for in F r (q). Three ranges of characteristic scales are 
observed. In the large-scale regime of r I, the mix- 
ture appears to be homogeneous. This leads that Eq.([3]) 
holds, and then F q (r) — 1 is obtained irrespective of 
q. At the small-scale regime of r <C I, the mixture has 
a certain level of inhomogeneity. In this scale, Eq.© 
holds and F q {r) - {A/A c ') 1/{q ' 1) = A' 1 = est. is ob- 
tained asymptotically irrespective of q. Thus, F q (r) ap- 
proaches a certain level depending on non-uniformity at 
the smallest scale. At the correlation scale of r ~ I, 
the rapid variation of F q (r) characterizes the correlation 
scale in the mixture. From the inhomogeneity function, 
we can identify both the macroscopic homogeneity and 
microscopic inhomogeneity of a mixture. To characterize 
macroscopic homogeneity, we define the scale of homo- 
geneity, Ih-, which is a lower bound, as F q (r) w 1 holds 
for Ih < r. 

In general, F q (r) is a monotonously decreasing function 
of the scale r, which corresponds to the decreasing inho- 
mogeneity as the observation scale enlarges. To identify 
a characteristic scale in a mixture clearly, it is convenient 
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to see the first derivative of F q (r). We define the scale 
susceptibility of the concentration fluctuation as 



G q {r) ee - 



_ dF q (r) 



d In i 



(8) 



which takes a finite positive value when F q (r) decreases 
sharply based on the characteristic scale of the fluctua- 
tion. Because we are foucusing on hierarchical structures 
with a wide range of scales, the scale susceptibility is 
defined as the derivative with respect to lnr. 

For the cases described in Eq. (0 , the scale susceptibil- 
ity becomes 



r < I, 
G q {r) — > <^ finite r ~ I, 
r > I, 



(9) 



from which the characteristic scale I is identified. In prac- 
tical situations, the scale susceptibility G q (r) is more con- 
venient than F q (r) to identify the characteristic scale I. 

We demonstrate how F q (r) and G q (r) work by apply- 
ing them to synthetic images with a definite characteris- 
tic scale and different levels of concentration fluctuation. 
We consider the model systems in Fig. [TJ Each system 
in Fig. Q] consists of 168x168 pixels, which define the 
range of r in the unit of pixel is from 1 to 168. and 
the concentration of a component is represented in gray- 
scale where the value of concentration, c(x), is within 
[0,255]. The single characteristic scale I of each system 
is naturally defined as the linear sizes of the darker ar- 
eas, which are 84, 28, 7 and 2 pixels from the left to the 
right columns. The systems in the first row model the 
strongly segregated mixtures where the dichotomous val- 
ues of concentration are 255 and 0. By moving down the 
matrix, the concentration fluctuation becomes weaker as 
the contrast between "black" and "white" reduces. 

The inhomogeneity functions of q = 2 for the systems 
in Fig. Q] are drawn in Fig. [5] Because the concentra- 
tion value is dichotomous in these systems, a single value 
of q is sufficient. Different values of q did not change 
the qualitative aspect of F q {r). The F q (r) for all the 
systems shows the sharp variation in each characteristic 
scale. The differences in the concentration fluctuation 
are explained by the magnitude of F q (r) at r <C I; the 
stronger the concentration fluctuation, the larger the in- 
homogeneity function. For this dichotomous- valued field, 
an F q {r — y 1) is at the finest resolution and is expressed as 

W) = 2 [{1 - (#)} g + (#)f , where W de- 
notes the intensity of "white" region. This result shows 
that the intensity of segregation is characterized by the 
value of F q (r) > 1 . 

The characteristic scale and intensity of segregation are 
independent characteristics of a mixture and are simul- 
taneously accounted for by F q {r). To isolate the char- 
acteristic scales, the scale susceptibility Gzir) for pat- 
terns in Fig. [T] is shown in Fig. [3] The G2(r) for the 
same scale of segregation is almost identical, except for 



the ordinate scale, and clearly indicates I. The useful- 
ness of G%{r) is for systems with small fluctuations but 
a definite geometric structure, such as the systems of 
("white", "black" ) = (127:128) in Fig.[Q The inhomogene- 
ity functions in Fig. [2fe) for these (127:128)-systems are 
close to unity at all scales, and the characteristic scales 
are almost indiscernible without magnification. However, 
the peaks in G2(r) in Fig. [3X e ) clearly show the I even 
for this weakly segregated system. 

We add some comments on the relationship between 
mixture quantification without accounting for a scale and 
the method of the scale-dependent moment. In the limit 
of r — ¥ 0, the moment function of a coarse-grained con- 
centration becomes (/^->o) — ^ In ^xp, g (x), which is pro- 
portional to the usual moment of the concentration field. 
The variance-based measures, such as the intensity of seg- 
regation, are basically related to this limit of q = 2. Such 
quantifications describe some average deviation from the 
uniform distribution in the smallest resolution and there- 
fore are insensitive to the scale of any internal structure. 




FIG. 1. Model systems of partially miscible blends of 168 x 168 
pixels: checkerboard patterns with varying shades of gray and 
varying characteristic scales of 84x84, 28x28, 12x12, 7x7, 
and 2x2 pixels. The intensities of black and white regions are 
the same in each row and are indicated in the first column. 



III. RESULTS AND DISCUSSION 

A. Application to particle dispersions 

A particle dispersion and a composite material are typ- 
ical examples of a multi-component system, which con- 
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FIG. 2. Inhomogeneity functions for the patterns in Fig. [T] as 
a function of the resolution scale. The "black" and "white" 
intensities are (a) 0:255, (b) 31:224, (c) 63:192, (d) 95:160 and 
(e) 127:128. 



sists of small particles and a matrix medium. The phys- 
ical properties of a dispersion are not solely determined 
by those of the particles and the matrix medium but 
also depend on the distribution of the particles. The 
structure of the particle distribution is controlled by the 
non-equilibrium processing history, as well as by the ther- 
modynamic stability of the structure. The mixture state 
of a dispersion should not be characterized solely by the 
primary particle size but also by multiple scales associ- 
ated with particle distributions. We apply the moment 
function method for characterization of different states 
of synthetic particle dispersions. 

Consider monodisperse dispersions with a diameter of 
the particles d = 10 pixels that is a natural character- 
istic scale of the systems. For the sake of simplicity, 
two-dimensional systems are analyzed, but application 
to three-dimensional systems is straightforward. Three 
different dispersions with an area fraction of the parti- 
cles of 0.1 are depicted in Fig. Hta)-(c). The particles 
are randomly distributed in Fig. HJa) , the particles are 
locally ordered in Fig. 0|b) and several particles aggre- 




r (pixel) 
(c) (White, Black)=(63, 192) 



r (pixel) 
(d) (White,Black)=(95,160) 




FIG. 3. Scale susceptibility to the resolution scale, G<i(r), for 
the systems in Fig.[T]as a function of the resolution scale. The 
"black" and "white" intensities are (a) 0:255, (b) 31:224, (c) 
63:192, (d) 95:160 and (e) 127:128. 



gate in Fig. HJc). The inhomogeneity function, F^ir), 
and the scale susceptibility, G&(r), for the three systems 
in Figs. |U[a)-(c) are shown in Fig. and Fig. BJe), 
respectively. In this case, the concentration field of the 
particles was analyzed, which was unity in the particle 
domain and zero in the matrix domain. Because the anal- 
ysis was applied to binarized systems, a single value of q is 
sufficient. For clarity of presentation, the results of q = 4 
are shown. For the random dispersion in Fig. f^a), there 
are no coherence scales other than d because the parti- 
cles are randomly distributed by construction. The F4 (r) 
and G±{r) for the random dispersion indicate the parti- 
cle size by the large fluctuation at r — d and the scale of 
homogeneity for £™ nd ° m « lfjd. The ratio l^ ndom /d = 10 
is relatively large, which explains the inhomogeneity of 
the particle distribution. 

In the system of Fig. 0|b), inter-particle distances be- 
tween nearest pairs are almost the same, and the parti- 
cles are locally ordered; therefore we call the system an 
ordered dispersion for convenience. The average inter- 
particle distance is r\ « 2.7 'd, which was measured by the 
first peak in the correlation function or the pair correla- 
tion function. The r\ and d are the natural characteristic 
scales in the ordered dispersion. The large fluctuation at 
r = d in F±(r) and G±{r) for the ordered dispersion indi- 
cates the particle size. At r > r\ in F±(r), the fluctuation 
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FIG. 4. Particle dispersion and their multiple-scale inhomo- 
geneity. Three different particle dispersions of 512x512 pix- 
els: (a) particles are randomly distributed, (b) particles are 
locally ordered, and (c) particles are aggregated. The area 
fraction of the systems is 0.1 and the diameter of the parti- 
cles is d — 10 pixels, (d) The inhomogeneity functions, F^r), 
and the scale susceptibilities, Ga(t), for the systems (a)-(c). 
Coarse-grained fields, H22(x), for (f) random dispersion, (g) 
ordered dispersion and (h) aggregated dispersion. 



almost vanishes because the structural fluctuation does 
not exists at large scale. A close examination of Gi(r) 
reveals a small peak at around r w 4<i, which can be 
ascribed to the fluctuation in a middle-range scale. The 
scale of homogeneity is about /™ dorod rj Qd. If the ordered 
configuration developed globally, ^ dcrcd « n would be 
expected because the fluctuation r > r\ did not exist. 
The ordered configuration makes the scale of homogene- 
ity smaller than that of the random dispersion. 

In the system of Fig. Etc), several particles aggregate, 
and these aggregates are distributed randomly in an ag- 
gregated dispersion. One aggregate is composed of 2-20 
particles, thus the aggregated dispersion should be char- 
acterized by the multiple sizes of the aggregates and the 
nearest- neighbor distance r\ = 1.2d in an aggregate be- 
sides the particle size d. In addition to the large fluctua- 
tion at r — d by the particle size, another broad band at 
n < r < lOd in G^{r) are observed, which corresponds to 
the multiple sizes of the aggregates. The scale of homo- 



geneity is about ^ grcgatc ~ ^Od, which is similar to that 
in the random dispersion and much greater than that 
in the ordered dispersion because of the random distri- 
bution of the aggregates. Difference of mixture quality 
among three dispersions is obvious in the scale suscep- 
tibility function, G^(r), in Fig. SJe). Identification of 
multiple characteristic scales based on inhomogeneity is 
essential to quantify mixture quality. 

Figures |4]T)-(h) show the coarse-grained fields of the 
particle dispersions in Figs. Hfa)-(c) with the resolution 
scale of r = 2. 2d. At this resolution scale, a single parti- 
cle is hardly discernible. In Fig. HJb), the fi22(x) for the 
ordered dispersion, which has the smallest Ih, shows the 
smallest contrast indicating the smallest fluctuation. For 
the random dispersion, fJ,22(x) shows a random structure 
in the fluctuation (Fig. 0Ja)). For the aggregated dis- 
persion, (J-22(x) shows a streak structure and the largest 
fluctuation among the three systems (Fig. Hfc)). 

Coarse-graining is a key concept in simultaneous char- 
acterization of various characteristic scales associated 
with a hierarchical structure and inhomogeneity of distri- 
bution. The demonstration with the different dispersions 
above showed how the scale-dependent moment method 
works to assess different characteristic scales including 
the particle size, aggregate sizes, and inter-particle spac- 
ing. 



IV. CONCLUDING REMARKS 

We proposed a method for characterization of hierar- 
chical mixture quality of complex mixtures , which we 
call the scale-dependent moment function method. In 
this method, the inhomogeneity of a mixture is defined 
at different scales. From this scale-dependent inhomo- 
geneity, multiple characteristic scales in mixture quality 
are identified. In other words, the sizes and distribu- 
tions of internal structures are simultaneously accounted 
for in the scale-dependent moment function. Application 
of this method to dispersions with different structures 
showed that differences in mixture quality are charac- 
terized by multiple length scales and associated inhomo- 
geneity. 

The scale-dependent moment method is applicable to 
two- or three- dimensional concentration fields. The mix- 
ture quality characterization based on the scales and 
inhomogeneity should be effective not only for multi- 
component complex mixtures with hierarchical internal 
structures but also for the evolution of the mixing pro- 
cess of both miscible and immiscible mixtures. It is ex- 
pected that this proposed multiple-scale characterization 
can provide an insight on the relationship between the 
distribution of internal structures and the physical prop- 
erties of mixture materials, but this issue is left for future 
work. 
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